Generating random density matrices 
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(Dated: March 18, 2011) 

We study various methods to generate ensembles of random density matrices of a fixed size N, 
obtained by partial trace of pure states on composite systems. Structured ensembles of random pure 
states, invariant with respect to local unitary transformations are introduced. To analyze statistical 
properties of quantum entanglement in bi-partite systems we analyze the distribution of Schmidt 
coefficients of random pure states. Such a distribution is derived in the case of a superposition 
of k random maximally entangled states. For another ensemble, obtained by performing selective 
measurements in a maximally entangled basis on a multi-partite system, we show that this distri- 
bution is given by the Fuss-Catalan law and find the average entanglement entropy. A more general 
class of structured ensembles proposed, containing also the case of Bures, forms an extension of the 
standard ensemble of structureless random pure states, described asymptotically, as N — > oo, by 
the Marchenko-Pastur distribution. 



I. INTRODUCTION 



Random states are often used in various problems in 
quantum mechanics and the theory of quantum informa- 
tion. On one hand, to describe properties of a quantum 
state affected by noise present in the system, one may 
assume that a given state is subjected to random interac- 
tion. On the other hand, random quantum states emerge 
in a natural way due to time evolution of arbitrary initial 
states of quantum analogues of classically chaotic systems 
[l[ . Furthermore, not knowing much about a given physi- 
cal state one can ask about their generic properties, char- 
acteristic of a 'typical' state. For instance, a key conjec- 
ture in the theory of quantum information concerning the 
additivity of minimal output entropy was recently shown 
to be false by investigating properties of random states 
obtained by random operations applied to the maximally 
entangled state of a composed system 0] . 

A standard ensemble of random pure states of size N 
is induced by the Haar measure over the unitary group 
U(N). The same construction works for quantum com- 
posite systems. For instance, random pure states of 
an N x K quantum system corresponds to the natural 
Fubini-Study measure, invariant with respect to the uni- 
tary group U (NK). Thus such ensembles of random pure 
states are structureless, as the probability measure is de- 
termined by the total dimension of the Hilbert space and 
is does not depend on the tensor product structure 0, 0] . 
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In this work we are going to analyze structured ensem- 
bles of random states on a composite systems, in which 
such a tensor product structure plays a crucial role. For 
instance, in the case of an N x K system, it is natural 
to consider ensembles of random states invariant with re- 
spect to local unitary transformations, described by the 
product unitary group, U(N) x U(K). Well known con- 
structions of random product states and random maxi- 
mally entangled states can thus serve as simplest exam- 
ples of the structured ensembles of random states. Other 
examples of structured ensembles of random pure states, 
which correspond to certain graphs were recently studied 
in@. 

The main aim of this paper is to introduce physically 
motivated ensembles of structured random states and to 
characterize the quantum entanglement of such states. 
For this purpose we analyze the spectral density P(x) of 
the reduced density matrix of size N, where x denotes 
the rescaled eigenvalue, x = NX. In some cases we eval- 
uate also the average entanglement entropy, defined as 
the Shannon entropy of the vector A of the Schmidt co- 
efficients. 

We treat in detail two cases of a direct importance 
in the theory of quantum information. The first ensem- 
ble is obtained by taking a coherent superposition of a 
given number of k independent, random maximally en- 
tangled states of an A x A system. We derive explicit 
formulae for the asymptotic spectral density Pk (x) of the 
reduced state and show that in the limit of large k, the 
density converges to the Marchenko-Pastur distribution 
(MP) [1, 0]- This feature is characteristic to the case of 
the structureless ensemble of random pure states, which 
incihil^esip^bato^lifcyllii<sfiujat€S\ta4ii.e set of mixed quantum 
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states dHH]. 

To introduce the second example of a structured en- 
semble we need to consider a four partite system. We 
start with an arbitrary product state and by allowing 
for a generic bi-partite interaction we create two ran- 
dom states on subsystems AB and CD. The key step 
is now to perform an orthogonal selective measurement 
in the maximally entangled basis on subsystems B and 
C. Even though the resulting pure state on the remain- 
ing subsystems A and D does depend on the outcome of 
the measurement, its statistical properties do not, and 
we show that the corresponding level density is given by 
the Fuss-Catalan distribution. 

Note that the above two constructions could be exper- 
imentally accessible, at least in the two-qubit case for 
k = 2. The second construction can be easily generalized 
for a system consisting of 2s subsystems. An initially 
product state is then transformed by a sequence of bi- 
partite interactions into a product of s bi-partite random 
states. Performing an orthogonal projection into a prod- 
uct of (s— 1) maximally entangled bases we arrive with a 
random bipartite state of the structured ensemble defined 
in this way. The distribution of its Schmidt coefficients 
is shown to be asymptotically described by the Fuss- 
Catalan distribution of order s, since its moments are 
given by the generalized Fuss-Catalan numbers [lC)l - [l2j ]. 
This distribution can be considered as a generalization of 
the MP distribution which is obtained for s = 1. We con- 
clude this paper proposing a generalized, two-parameter 
ensemble of structured random states, which contains, as 
particular cases, all the ensembles analyzed earlier in this 
work. 

This paper is organized as follows. In section II we re- 
call necessary definitions and describe a general scheme 
of generating mixed states by taking random pure states 
on a composite system of a given ensemble and taking 
an average over selected subsystems. In Section III we 
use 2 independent random unitary operators to construct 
the arcsine ensemble obtained by superposition of two 
random maximally entangled states. More general struc- 
tured ensembles are obtained by superposing k random 
maximally entangled states. Other ensembles of random 
pure states are discussed in section IV. For completeness 
we review here some older results on ensembles which 
lead to the Hilbert-Schmidt and the Bures measures in 
the space of density matrices. 

In section V we present a scheme to generate random 
states by considering a state defined on a 2s-partite sys- 
tem and performing measurement onto the product of 
(s — 1) maximally entangled states. The distribution 
of the Schmidt coefficients of the resulting random pure 
state is then given by the Fuss-Catalan distribution of or- 
der s. An explicit form of this distribution is presented in 
Sec. VI for any s. In Sec. VII we list various ensembles 
of random states, characterize their statistical properties 
and introduce a generalized ensemble of structured states 
P2l . Details concerning the properties of the sum of 
k independent random unitary matrices are presented in 



Appendix A. Brief discussion of real random pure states 
and real density matrices is provided in Appendix B, in 
which we state that the distribution of singular values 
of a product of s real random Ginibre matrices can be 
asymptotically described by the Fuss-Catalan distribu- 
tion. 



II. ENSEMBLES OF RANDOM STATES 

To construct an ensemble of random states of a given 
size N one needs to specify a probability measure in the 
set of all density operators of this size. Interestingly, 
there is no single, distinguished probability measure in 
this set, so various ensembles of random density operators 
are used. 

A possible way to define such an ensemble is to take 
random pure states of a given ensemble on bi-partite 
system and to average over a chosen subsystem. The 
structureless ensemble of random pure states on N x N 
system distributed according to unitarily invariant mea- 
sure leads then to the Hilbert-Schmidt measure 0, 0]- A 
similar construction to generate random states according 
to the Bures measure Il4j was recently proposed 

in [l5j]. These algorithms to generate random quantum 
states are already implemented in a recent Mathematica 
package devoted to Quantum Information [l6|, [l?} ■ 

Let us recall here the necessary notions and defini- 
tions. Consider the set of quantum states M.n which 
contains all Hermitian, positive operators p = > 
of size N normalized by the trace condition Trp = f . 
Any hermitian matrix can be diagonalized, p — VAV* . 
Here V is a unitary matrix consisting of eigenvectors of p, 
while A is the diagonal matrix containing the eigenvalues 
(Ai, . . . , \n}- 

In order to describe an ensemble of random density 
matrices, one needs to specify a concrete probability dis- 
tribution in the set A4n of quantum states. In this work 
we are going to analyze ensembles of random states, for 
which the probability measure has a product form and 
may be factorized 0, [H| , 

dp, p = di/(Ai, A 2 ,..., Ajv) X d/zy, (1) 

so the distribution of eigenvalues and eigenvectors are in- 
dependent. It is natural to assume that the eigenvectors 
are distributed according to the unique, unitarily invari- 
ant, Haar measure dpv on U (N). Taking this assumption 
as granted, the measure in the space of density matrices 
will be determined by the first factor dv describing the 
joint distribution of eigenvalues P(Ai, . . . , Ajv). 

In quantum theory mixed states arise due to an interac- 
tion of the system investigated with an external environ- 
ment. One may then make certain assumptions concern- 
ing the distribution of pure states describing the extended 
system. The desired mixed state p on the principal sys- 
tem A of size N can be then obtained by partial trace 
over the subsystem B of an arbitrary size K. 
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Consider an arbitrary orthonormal product basis \i) ® 
\j) E H N ®H K , with i = 1, .. . , N and j = 1,...,K. Any 
pure state \ip) of the bi-partite system (not necessarily 
normalized) can be expanded in this base, 

N K 

w = ®ii>> ( 2 ) 

i=i j=i 

where X is a given complex rectangular matrix of size 
JV x if. 

Let us then consider an arbitrary complex matrix X. 
It leads to a (weakly) positive matrix XXK Thus nor- 
malizing it one obtains a legitimate quantum state, which 
corresponds to the partial trace of the initial pure state 
\t/j) over the subsystem B, 

P (V#) TrAAt ■ [ ' 

For instance, taking X from the Ginibre ensemble fl8l.ll9l] 
of complex square matrices we get the Hilbert-Schmidt 
ensemble of quantum states, while a more general family 
of induced measures corresponds 20] to the ensemble of 
rectangular Ginibre matrices. 

The spectrum of a density matrix p is thus equivalent 
to the set of Schmidt coefficient of the initially pure state 
\ip), which are equal to squared singular values of a ma- 
trix X, normalized in such a way that Trp = 1. The 
degree of mixing of the reduced matrix p can be charac- 
terized by its von Neumann entropy, S(p) = — Trp In p, 
equal to the Shannon entropy of the Schmidt vector, 
-^'(IV')) = ~~^2i^i^ n ^i- This quantity is also called en- 
tropy of entanglement of the pure state \ip), as it is equal 
to zero iff the state has a tensor product structure and 
is separable. We are also going to use the Chebyshev 
entropy, which depends only on the largest eigenvalue 
Soa{p) = — InAmax, and determines the geometric mea- 
sure of entanglement of the pure state \ip). 

Different assumptions concerning the distributions of 
pure states of the bi-partite system lead to different en- 
sembles of mixed states on the system A. A natural as- 
sumption that \tp) belongs to the standard, structureless 
ensemble of random pure state distributed with respect 
to the unitarily invariant measure leads to the induced 
measures [4j in the space of mixed states. 

In the subsequent section we are going to present sim- 
ple examples of structured ensembles of random states, 
in which the tensor structure plays a key role. For in- 
stance, we discuss first ensembles of random states of a 
bi-partite systems, which are invariant with respect to 
local unitary operations. 



III. TWO PARTITE SYSTEMS 

A. Random separable and maximally entangled 
pure states 

For completeness we shall start the discussion with 
the somewhat trivial case of generating random sepa- 
rable states. Consider an arbitrary product state on a 
bi-partite system, |0, 0) = |0)a ® |0)b- A local unitary 
operation, U — Ua ® Ub, cannot produce quantum en- 
tanglement, so the state defined by two random unitary 
matrices, \iPab) = U\0,0) = Ua\0)a®Ub\0)b is also sep- 
arable. Hence a random separable state is just a product 
of two random states, \i/>ab) = \iPa) ® IV'fl) 

Consider now a generalized Bell state defined on an 
N x N system A, B, 

1 N 
v i=l 

As the entanglement entropy of this state is maximal, 
E(\^> + )) = In N, this state is called maximally entangled. 

Any local unitary operation, Ua <8> Ub, preserves 
quantum entanglement, so the locally transformed state 
\4>cnt) '■= (Ua ® Ub)\^ + ) remains maximally entangled. 
Hence taking random unitary matrices Ua,Ub £ U(N) 
according to the Haar measure we obtain an ensemble of 
random entangled states, such that their partial trace is 
equal to the maximally mixed state, 

p = Tr B |0 cnt )(0 cnt | = -t^Ijv =: P*- (5) 

This Dirac distribution will be denoted by 7r(°) (p) = 5(p— 
p*) and a scheme of generating it by averaging over an 
auxiliary subsystem is shown in Fig. [T] 

Pa Tr B 

FIG. 1: Generating mixed states according to ty^: step i) 
a local operation Ua <S> Ub on |^^ B ) creates a random en- 
tangled state |<^cnt), while the partial trace over an auxiliary 
subsystem B leads in step ii) to the maximally mixed state 
on the system A. 



B. Arcsine ensemble 

Consider now a more general case of a coherent super- 
position of two maximally entangled states. To be precise 
we fix a given maximally entangled state IV'i), and use 
such a local basis that it is represented by (j4|). Taking 
a local random unitary matrix Ua € U(N) we generate 
another maximally entangled state = Ua ® 
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As shown in Fig. [2] we construct their symmetric super- 
position, 

\<t>) = + life}) = [liv= +{U® 1nW + ) ■ (6) 

Note that this ensemble is invariant with respect to lo- 
cal unitary operations, U(N) x U(N). Let us specify a 
subsystem B and average over it to obtain the reduced 
state 



P = 



21 + U + W 
2N + Tr(U + W)' 



(7) 



If the matrix U is generated according to the Haar 
measure on U (N) its eigenphases are distributed accord- 
ing to the uniform distribution 19], P(a) — l/2n for 
a e [0,2tt). Thus the term Tr(U + W) present in the 
normalization constant becomes negligible for large N . 
Therefore eigenvalues of a random density matrix p by 
(0 have the form Aj = (1 + cos ai)/N, for i = 1, . . . , N. 
Making use of the rescaled variable x = XN we arrive 
at a conclusion that the spectral level density of random 
density matrices ([7]) is asymptotically described by the 
arcsine distribution 21] 



Pare 0*0 



1 



TT^/x(2 



(8) 



defined on the compact support [0,2]. This is a par- 
ticular case of the beta distribution (with parameters 
a = /3 = 1/2) and the Dirichlet distribution. It de- 
scribes the Jeffreys prior of a Bernoulli trial [22j and is 
related to the statistical distance between classical prob- 
ability distributions Q. The name 'arcsine' is due to 
the fact that the corresponding cumulative distribution 
is proportional to sin -1 yJx/2. Thus the ensemble of 
random density matrices constructed according to the 
procedure shown in Fig. [5] will be called arcsine ensem- 
ble. The average entropy for the arcsine distribution 
reads / Q 2 — x\nxP a rc{x) = In 2 — 1 « —0.307, so the aver- 
age entropy of entanglement of a random pure state on 
the N x N system generated form this ensemble reads 



InN- 1 + ln2. 



l*AB + > 
• <r > 



B 




u 



A 



FIG. 2: To generate states from the arcsine ensemble de- 
scribed by (0 one has to i) construct a superposition of a 
maximally entangled state |<&^ s ) with another maximally en- 
tangled state (Ua®1)\ } $\ b } , and ii) perform partial trace over 
an auxiliary subsystem B. 

It is clear that one may define a family of interpolating 
ensembles by taking a convex combination of the mixed 



states defined by (J5J) and 0. In other words, one can 
take a family of random matrices parametrized by a real 
number a € [0, 1] and write W a = ol+(l— a)U. Plugging 
this expression in place of X into © we construct a fam- 
ily of ensembles of density matrices which gives the Dirac 
mass for a = and a = 1, while the arcsine ensemble is 
obtained for a — 1/2. 



C. Superposition of k maximally entangled states 

The arcsine ensemble introduced above can be ob- 
tained by superimposing k = 2 random maximally en- 
tangled states. It is straightforward to generalize this 
construction for an arbitrary number of k maximally en- 
tangled states. Each of them can be written by an action 
of a local unitary on the fixed maximally entangled state 
|*+). More precisely we set \ipi) = (Uj ® 1)|* + ) for 
i = 1 , . . . , k and construct their equi-probable superposi- 
tion |0) = Eti = [(Eli Ui) ® In] |*+>. 

As before this ensemble is invariant with respect to 
local transformations. The random mixed state is ob- 
tained by taking the partial trace over the subsystem B 
and normalizing the outcome p — Ti b\4>) (4>\/ (4>\4>) as in 
cq. ([7]). This procedure leads now to the following ex- 
pression for a mixed state, generated by k independent 
random unitaries, 



u k )(u} 



ul) 



u k ){ul 



(9) 



As shown in Appendix A the spectral density of ran- 
dom states defined above converges, for large system 
size N to the following distribution Vk{x) defined for 



Vk(x) 



1 y/Ak(k- l)x-k 2 x 2 
2n kx — x 2 



(10) 



The shape of these distributions is presented for some 
values of k in Fig. |3l It is easy to see that for a large num- 
ber k the above measures converge weakly to the limit 
Vac = = ^ A x-x 2 which i s the Marchenko-Pastur dis- 
tribution tt^. In other words, the superposition of a 
large number of random maximally entangled states de- 
stroys the structure of the ensemble, as the resulting state 
becomes typical of the structureless ensemble. 

The measure (|10[) can be connected to the free Meixner 
measures of [23|: v k is the free Meixner measure of pa- 
rameters a — and b = —1/k. The measure v k {x) 
can be characterized by its second moment, which gives 
the asymptotic average purity of random mixed states 
obtained by a superposition of k maximally entangled 
states, (Trp 2 )^ w (2- l/fc)^. 




FIG. 3: Family of spectral distributions Vk(x) plotted for 
k = 2, 3, 5, 10 and k — > oo. The support of the measures is an 
increasing function of the parameter k. They describe level 
density of mixed states generated by superposition of k ran- 
dom maximally entangled states and performing the partial 
trace. The case k = 2 corresponds to the arcsine ensemble, 
while in the case k — > oo the sum of k random unitaries has 
properties of a random Ginibre matrix, so the spectral distri- 
bution converges to the Marchenko-Pastur law. 



D. Induced measures 

Consider a quantum system composed of two subsys- 
tems. Let N denote the dimension of the principal sys- 
tem, which interacts with a K dimensional environment. 
Assume first that the state £ Hn^Hk is taken from 
the structureless ensemble, so it is distributed uniformly 
according to the Fubini-Study measure on the complex 
projective space. In other words, the random pure state 
can be represented by \tp) = U\Q, 0), where U is a global 
random unitary matrix distributed according to the Haar 
measure - see Fig. 2J The initial state is arbitrary, so for 
concreteness we may choose it as a given product state 
|0, 0) in Hn <8> Hk- The mixed state p obtained by the 
partial trace over the K dimensional environment 



p = Thfl^X^I , with \t/j)eH N <S>H 



(11) 



is distributed according to the measure pn.k in the space 
of density matrices of size N, induced by the Haar mea- 
sure on the unitary group U(NK). This procedure yields 
the following probability distribution 



dfx(p) cx @(p) 5{Trp - 1) det p 



K-N 



(12) 



The O step-function and the Dirac S reflect key proper- 
ties of density matrices: positivity, p > and normaliza- 
tion, Trp = 1. 

If is a random pure state distributed uniformly 
according to the unitarily invariant Fubini-Study mea- 
sure in the space of pure states all elements of G are 




'A 



Tr 



13 



11) 




FIG. 4: Taking i) a random pure state Uab\0, 0} on a bi- 
partite system AB and performing ii) the partial trace over 
subsystem B leads random state distributed according to 
an induced measure, with spectral density described by the 
Marchenko-Pastur distribution. 



independent complex Gaussian variables with the same 
variance, so this matrix belongs to the Ginibre ensemble 
[l8j |. Hence the reduced state (|TT|) obtained by partial 
trace has the form with X = G. The eigenvalues of a 
random matrix p generated with respect to the measure 
Pm,k are thus equal to the squared singular eigenvalues 
of a normalized rectangular N x K complex random ma- 
trix G from the Ginibre ensemble. 

Any density matrix is Hcrmitian and can be diagonal- 
ized by a unitary rotation. Integrating out the eigenvec- 
tors of a random state p defined by (fT2|) one reduces dp, to 
the measure on the simplex of eigenvalues {Ai, . . . , \n} 
of the density operator [24j , 

1...JV 

Pn,k(\u-\n) = C n , k Y[\?- n HiK-Xj) 2 . (13) 

% i<j 

In this expression all eigenvalues are assumed to be non- 
negative, Ai > and they sum to unity, J^Aj = 1. 
The normalization constant can be expressed [4[ by the 
Gamma function 



Cn,k — 



T(KN) 



■j)T(N-j + l) 



(14) 



Induced measures are stable with respect to the partial 
trace, what can be formulated as follows. 

Proposition 1. Consider a random state p gener- 
ated according to the induced measures pm,k- Assume 
that the dimension is composed, M — N x L, and de- 
fine the partial trace over the L dimensional subsystem, 
a = Ttl(p). Then the reduced matrix a is generated 
according to the measure pn,kl- 

Proof. The state p can be considered as a random 
pure state \ip) € W.m®T-Ik averaged over K dimensional 
environment. Averaging the projector IV'XV'I over a com- 
posed subsystem of size KL we arrive at state a of size 
N, which shows [25[ that it is distributed according to 
the induced measure pn,kl- 



IV. MEASURES DEFINED BY A METRIC 

Let d denote any distance in the space of the normal- 
ized density operators of a fixed size N. With respect to 
this distance one can define unit spheres and unit balls. 
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Attributing the same weight to any ball of radius equal 
to unity one defines the measure corresponding to a given 
distance. 

Consider, for instance, the Hilbert-Schmidt distance 
between any two mixed states Dhs{p,&) = [Tr(p — 
ct) 2 ] 1 / 2 . This distance induces in the space of density op- 
erators an Euclidean geometry: the space of one-qubit 
mixed states has a form of the three-ball bounded by the 
Bloch sphere. 



As the rescaled eigenvalue is x = NX, this result is consis- 
tent with the fact that the average entropy of the asymp- 
totic distribution (ITTl) reads J Q 4 — xlnxPns(x)dx = —1/2. 
Note that the mean entropy of a random mixed state is 
close to the maximal entropy in the N dimensional sys- 
tem, which equals IniV for the maximally mixed state 
= 1/N. The asymptotic average purity (Trp 2 )HS ~ 
2/N which is consistent with the second moment of the 
MP distribution, J*x 2 Prs(x) = 2. 



Hilbert-Schmidt ensemble 



B. Bures ensemble 



The Hilbert-Schmidt measure denned in this way be- 
longs to the class of induced measures [3| and can be ob- 
tained by a reduction of random pure states denned on a 
bi-partite quantum system. Looking at expression (|13[) 
we see that in the special case K = N the term with the 
determinant in is equal to unity, so measure reduces to 
the Hilbert-Schmidt measure |J|. This observation leads 
to a simple algorithm to generate a Hilbert-Schmidt ran- 
dom matrix [4J : i) Take a square complex random matrix 
A of size N pertaining to the Ginibre ensemble [HI, [l9[ 
(with real and imaginary parts of each element being in- 
dependent normal random variables); ii) Write down the 
random matrix 



Phs 



TrGGt ' 



(15) 



which is by construction Hermitian, positive definite and 
normalized, so it forms a legitimate density matrix. 

To characterize spectral properties of random density 
matrices one considers the joint distribution of eigenval- 
ues (fT5|) , integrates out N — 1 variables Ai . . . Ajv to ob- 
tain the level density denoted by P(X). This problem 
was discussed by Page [|[ , who found the asymptotic dis- 
tribution for the rescaled variable x — NX. The result 
depends on the ratio of both dimensions c = K/N and is 
given by the Marchenko-Pastur distribution [6] , 



(x) = max(l - c, 0)5(0) + 



1) 



2nx 



(16) 



This expression is valid for x £ [x—,x+], where x± = 
1 + c ± 2-y/c. In the case c = 1 corresponding to the 
Hilbert-Schmidt ensemble this distribution reads 



Pns(x) = ^J--l for ze[0,4] (17) 

27T V X 

It diverges as x -1 / 2 for x — > and in the rescaled vari- 
able, y — \fx it becomes a quarter-circle law, P(y) = 

The average von Neumann entropy of random states 
distributed with respect to the HS measure reads @, H(| 



1 /In N 



(18) 



Another distinguished measure in the space fl of quan- 
tum mixed states, is induced by the Bures distance 

mm, 



D B {p,a) = J 2 - 2Tr(VWp) 



1/2 



(19) 



The Bures distance plays an important role in the investi- 
gation of the set of quantum states The Bures metric, 
related to quantum distinguishability 22] , is known to be 
the minimal monotone metric |28| and applied to any two 
diagonal matrices it gives their statistical distance. These 
special features support the claim that without any prior 
knowledge on a certain state acting on Hn, the optimal 
way to mimic it is to generate a random density operator 
with respect to the Bures measure. 

The Bures measure is characterized by the following 
joint probability of eigenvalues [l3[ 

P B (X 1 ,...X N ) = C§l[xr^ J] %=^-. (20) 



i<3 



Xi + Xj 



where all eigenvalues are non-negative, Ai > and they 
sum to unity, J^. A,- = 1. The normalization constant for 
this measure 



C B 9 N*-N T(iV 2 /2) 



n N/2 n; =1 r( J + i) 



(21) 



was obtained in [3, [lij for small N and in [2{| in the 
general case. 

To generate random states with respect to the Bures 
ensemble one can proceed according to the following al- 
gorithm [rsl j: 

i) Take a complex random matrix G of size N pertain- 
ing to the Ginibre ensemble and a random unitary matrix 
U distributed according to the Haar measure on U(N) 

m. 

ii) Write down the random matrix 



PB 



(1 + U)GG\1 + W) 
Tr[(l + [/)GGt(i + [/t)] 



(22) 



which is distributed according to the Bures measure. 
In the analogy to the Hilbert-Schmidt ensemble we can 
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write this state as reduction of a pure state on the com- 
posed system, 



PROJECTION ONTO MAXIMALLY 
ENTANGLED STATES 



PB = where |0>:=[(1 + Va)®1]|V>i>. 

(</#) 

(23) 

Here 1^1)= C^ab|0,0) is a random state on the bipartite 
system used in Eq. {TT} and Va € U(N). 





FIG. 5: To generate states according to the Bures distribution 
one i) constructs a superposition of a random bi-partite state 
\ipi) = Uab\0,0) with a locally transformed random state 
~ (Va £S> l)\tpi), and ii) perform partial trace over an 
auxiliary subsystem B. 



The asymptotic probability distribution for the 
rescaled eigenvalue x — NX of a random density matrix 
generated according to the Bures ensemble reads [2Cj 




Pb(x) =C 



where C = 1/4W3 and a = 3\/3- This distribution is 
defined on a support larger than the standard MP distri- 
bution, x <E [0, a] and it diverges for x — » as £~ 2 / 3 . 

The average entropy of a random state form the Bures 
ensemble reads 120] 



(S(p)} B = In TV -hi 2 + 



In TV 



(25) 



This value is smaller than the average entropy (|18|) , which 
shows that the Bures states are typically less mixed that 
the states from the Hilbert-Schmidt ensemble. A similar 
conclusion follows from the comparison of average purity 
for the Bures ensemble, (Trp 2 ) B « 5/2iV ;20j, with the 
average purity for the HS measure. 

By considering random matrices of the form W = 
(at + (1 — a)U)G one may construct a continuous fam- 
ily of measures interpolating between the Bures and the 
Hilbert-Schmidt ensembles [15[ and labeled by a real pa- 
rameter a £ [0, 1/2]. A more general class of interpolat- 
ing ensembles is proposed in Sec. IVIII 



A. Four-partite systems and measurements in a 
maximally entangled basis 

Consider a system consisting of four subsystems, la- 
beled as A, B, C and D. For simplicity assume here that 
the dimensions of all subsystems are equal, N\ — N2 — 
N3 = N4 — N. Consider an arbitrary four-partite prod- 
uct state, saylV'o) = |0)A®|0) B <E>|0) c ®|0)rj =: |0, 0,0,0). 
Taking two independent random unitary matrices Uab 
and Ucd of size N 2 , which act on the first and the sec- 
ond pair of subsystems, respectively, we define a random 
state \ip) — Uab ® Ucd\4>o)- By construction, it is a 
product state with respect to the partition into two par- 
ties: (A,B) and (C,D). In the analogy to © it can be 
expanded in the product basis, 



N 



N 



W=J2 G V E u ® \3)b ® \ k )c ® \1)d (26) 

i,j=l k,l=l 

Consider now a maximally entangled state on the sec- 
ond and the third subsystem, 



\n c ) 



1 



N 

£lM>J 



® \n)t 



(27) 



and the corresponding projector Pbc '■= 



BC\ 



One can extend it into a four-partite operator and define 



P:= lx® (^53|m./*)bc<i/,H) ® 1 



(28) 



/1 , 1/ 



Let us assume that the random pure state \ip) defined in 
(1261) is subjected to a projective measurement performed 
onto the second and third subsystem, and that the result 
is post-selected to be associated to the projector P. This 
leads to a non- normalized pure state \<j>) describing the 
remaining two subsystems, 



N N N 



P\rl>) 



D ■ 



(29) 



i=i z=i fc=i 



\* + Bc) 



Note that the projection on the entangled state 
introduces a coupling between the subsystems B and C. 
A similar idea was used in analysis of entanglement swap- 
ping [34]]. and in studies of 'matrix product states' (3l| 
and 'projected entangled pair states' [32|,|33]. Construct- 
ing the matrix product states an N x N entangled state is 
projected down into a subspace of an arbitrary dimension 
d, while in our approach a projection onto the maximally 
entangled state of BC takes place, (which formally cor- 
responds to d = 1), and only two edge systems labeled 
by A and D survive the projection. 

Normalizing the resulting state \<fi) and performing the 
partial trace over the fourth subsystem D we arrive at a 
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compact expression for the resulting mixed state on the 
first subsystem A, 

T^dWM = GEE^G^ 
9 (<j>\<f,) Tr GE&Gft' { ' 

For any matrices G and E this expression provides a valid 
quantum state, normalized and positive. If initial pure 
states are random, then the matrices G and E belong 
to the Ginibre ensemble and the spectrum of a random 
state p consists of squared singular values of the product 
GE of two random Ginibre matrices. Random states 
with the same statistical properties were recently found 
in ensembles associated with certain graphs Q . 




BC. However, these results are equivalent up to a uni- 
tary transformation, which again can be absorbed into 
the definition of the Ginibre ensemble. Thus the statisti- 
cal properties of the random state ([50)1 on subsystem A 
obtained in consequence of the measurement in the maxi- 
mally entangled basis in BC followed by the partial trace 
over subsystem D do not dependent on the outcome of 
the measurement. 

Furthermore, the same construction holds in a more 
general setup in which some dimensions of four subsys- 
tems are different. To use the maximally entangled state 
\^bc) we set = Nc, but the dimensions N = Na 
and Nd can be different. This leads to two rectangular 
random Ginibre matrices, G of size N x Nb and E of 
size Nb x Nd. As in the previous case formula (f5U|) pro- 
vides a density matrix p of size N, but now the model is 
a function of two parameters: dimensions Nb and Nd- 
It is sometimes convenient to use two ratios, c\ — Nb/N 
and c 2 = Nd/N, so the standard version of the model 
corresponds to putting c\ = c 2 = 1. 



B. Multi— partite systems 



11) 



n 



A 



n 



B 



n 



c 



<- > 

l*BC + ><*BC + l 



n 



D 



Another possibility to generalize the model is to con- 
sider a larger system consisting of an even number 2s of 
subsystems. For simplicity assume first that their dimen- 
sions are set to N. In analogy to (j2"rJ)l we use s indepen- 
dent unitaries of size N 2 to generate a random pure state 



P 



A 



Tr 



D 



iii) 




FIG. 6: To generate a random mixed state with spectral 
density 7r' 2 ' one has to i) take a product random pure state 
\^>ab) ® \<Pcd) on a four particle system A, B,C, D, ii) mea- 
sure subsystems B, C by a projection onto the maximally en- 
tangled state iii) average over the subsystem D. 

Observe that the statistical properties of the ensem- 
ble defined will change if the projection is performed 
with respect to an arbitrary maximally entangled state, 
\^bcY = (.Ub®Uc)\^bc)> as the local unitaries Ub and 
Uc can be absorbed in the definition of random Ginibre 
matrices, G and E, respectively. 

Furthermore, one may chose N 2 unitary matrices Ui € 
U(N), which are orthogonal in sense of the Hilbert- 
Schmidt scalar product, Tr UW — NSij. Then the set of 

N 2 maximally entangled states, 1*+) = (Ui ® 1)|*bo)' 
forms a maximally entangled basis, which are known to 
exist in any dimension [35| . Thus one may consider an- 
other set up in which a selective measurement on subsys- 
tems B and C is performed in the maximally entangled 
basis. The outcome state on subsystems AD depends on 
the result of the quantum measurement of subsystems 



To work with an arbitrary even number of subsystems 
it is convenient to modify the notation and label the sub- 
systems by integers 1,2,..., 2s. Consider an arbitrary 
product state of a 2s - particle system, |^o) = |0)i ® 
• • • <E> |0)2s =: |0, . . . , 0). Taking s independent Haar ran- 
dom unitary matrices C/1,2, ^3,4, • • • C^2s-i,2s of size N 2 , 
we define a random state 1^) = C/1,2 ® • • • C/2 S -i,2s|^o}- 
Expanding this state in the product basis one obtains 

W - (C/1,2® •••C/ 2s -i, 2s )|0,...,0) (31) 



H,...ia 



(G 



1 )ii -Ai 



(G 



S/l2s-l)«2, 



■ ■ ■ Ms) 



Performing a projection onto a product of s— 1 maximally 
entangled states, 



Ps := ll 



WW®- 



)1 



2.s 



we obtain a pure state 
subsystems, 



(32) 

describing the remaining two 



= P\4>) =/V 1 - s ^(g 1 G 2 ---G s ) Ji)i®|i)a,. (33) 



Normalizing this state and performing the partial trace 
over the last subsystem we obtain an explicit expression 
for the resulting mixed state on the first subsystem 

Tr 2s |0)(0| G 1 G 2 ---G S (G 1 G 2 ---G s )i 

P = 



Tr [G 1 G 2 ---G S (G 1 G 2 ---G s )n 



(34) 



9 



Alternatively, one may assume that this state is obtained 
as a results of an orthogonal measurement into the prod- 
uct of s — 1 maximally entangled bases. The first entan- 
gled basis correlates subsystem 2 with subsystem 3, the 
next one couples subsystem 4 with 5, while due to the 
last one the subsystem 2s — 2 is correlated with 2s — 1. 
Thus eigenvalues of a random state generated in this way 
coincide with squared singular values of the product of s 
independent random Ginibre matrices. Their statistical 
properties will be analyzed in the following section. In 
general the Ginibre matrices need not to have the same 
dimension so the model can be generalized. Assuming 
that a rectangular matrix Gi has dimensions Ni x Mi 
one has to put Mj = Ni+i for i = 1, . . . , s— 1, so that the 
product l]34p is well defined. Setting Ni — N and defin- 
ing the ratios Cj = Mi/N for i = l,...,s one obtains 
ensemble of random states parametrized by the vector of 
coefficients, c := {ci, . . . c s }. 



i) 






U 



1,2 



u 



3,4 



TL\ Ti.2 7~^3 ^4 



ii) 




U 2s -l,2s 

^2s-l ^2s 
• • • 

<■-> 

P2s-2,2s-l 



lii) 



FIG. 7: To obtain a random mixed state with spectral density 
7T use a system consisting of 2s subsystems, 1, . . . , 2s of size 
N: i) take the product of s bi-partite random pure states gen- 
erated by random unitary matrices E/1,2, ^3,4, • • • U2 S -i,2s £ 
U(N 2 ), ii) measure subsystems 2, . . . , 2s — 1 by a projec- 



tion onto the product of s 

-P23 ® Pi,5 ® • ' ' ® P2s-2,2s- 



1 maximally entangled states, 



where 



perform partial trace average over the subsystem 2s. 



VI. PRODUCT OF GINIBRE MATRICES AND 
FUSS-CATALAN DISTRIBUTION 

For any integer number s, there exists a probability 
measure tt^ s \ called the Fuss-Catalan distribution of or- 
der s, whose moments are the generalized Fuss-Catalan 
numbers [T5L Il2j given in terms of the binomial symbol, 



b(s) 



X m TT {s) {x)d X 



1 / sm + m 
sm + 1 V m 



--: FC [s) 



(35) 

The measure 7r( s ) has no atoms, it is supported on 
[0, b(s)] where b(s) = (s + l) s+1 /s s , its density is analytic 



on (0, b(s)), and bounded at x = b(s), with asymptotic 
behavior ~ l/(nx s ^ s+1 ^) at x — > 0. This distribution 
arises in random matrix theory as one studies the prod- 
uct of s independent random square Ginibre matrices, 
W = Ylj=i Or I n tms case squared singular values of 
W (i.e. eigenvalues of W^W-^) have asymptotic distribu- 
tion tt( s ). The same Fuss-Catalan distribution describes 
asymptotically the statistics of singular values of s-th 
power of a single random Ginibre matrix [36| . In terms of 
free probability theory, it is the free multiplicative convo- 
lution product of s copies of the Marchenko-Pastur dis- 
tribution [ll|,[37j], which is written as tt^ = [ttW] 1 ^. 



0.25 




FIG. 8: Fuss-Catalan distributions ir^(x) plotted for s = 
1,2,3,4 supported on the interval [0, (s + l) s+1 /s s ]. To 
demonstrate the behavior at the right edge of the support 
the figure is depicted for x > 1. 

An explicit expression of the spectral density for s = 2, 



r(2) 



(X) 



S/2V3 [$2 (27 + 3V81 - Ux) 3 - 6$x\ 



12tt 



xi (27 + 3^81 - Ux) 



_ (36) 

where x £ [0, 27/4], was derived first in [38( in context of 
construction of generalized coherent states from combi- 
natorial sequences. More recently it was applied in [j| to 
describe random quantum states associated with certain 
graphs. 

The spectral distribution of a product of an arbitrary 
number of s random Ginibre matrices was recently an- 
alyzed by Burda et al. 39] also in the general case of 
rectangular matrices. The distribution was expressed as 
a result of a polynomial equation and it was conjectured 
that the finite size effects can be described by a simple 
m ultip licative correction. Another recent work of Liu et 
al. [40( provides an integral representation of the distribu- 
tion 7i"( s ) derived in the case of s square matrices of size 
N, which is assumed to be large. 

Making use of the inverse Mellin transform and the 
Meijer G-function one may find a more explicit form of 
this distribution. It can be represented 41( as a super- 
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position of s hypergeometric functions of the type s F s —x, 
I 



7r (a) (x) = ^ A„ jS x^r- 

n=l 

where the coefficients A„ 



A, 



read for n 



:= s- 3 / 2 



1,2,. ..,s 

a + 1 /s s /( s+1 K« 



2tt 



V s + 1 / 



rij=i ^(i+i) rij=n+i r( j s+1 ) 



(38) 



Here p-Fg ([{<%•} j = i] ; [{&j}j=i] > x ) stands for the hypergeometric function [42[ of the type p F q with p 'upper' pa- 
rameters cij and q 'lower' parameters bj of the argument x. The symbol {a,i} r i=1 represents the list of r elements, 
di, . . . a r - The above distribution is exact and it describes the density of squared singular values of s square Ginibre 
matrices in the limit of large matrix size N. 

Observe that in the simplest case s = 1 the above form reduces to the Marchenko-Pastur distribution, 



^ 1 Hx) = ^ 7 = 1 F ([-h[};l 

7TV X V 2 4 



TTVX 



while the case s = 2 



7r< a >(*) = 



V3 



2 '^ 

is equivalent to the form (|36|) obtained in [38j j . 



6' 3 J ' L 3 J ' 27 



V3 



^ / 1 2 4 4x\ 

V3 2 Fl V6' 3 J '^3 J; 27 J 



67ra; 



(39) 



(40) 



The distributions (|37]) are thus directly applicable to 
describe the level density of random mixed states ob- 
tained from a 2s-partite pure states by projection onto 
maximally entangled states and partial trace as described 
in previous section. This result becomes exact in the 
asymptotic limit, it the dimension N of a single subsys- 
tem tends to infinity. However, basing on recent results 
of Burda et al. [3S| one can conjecture that the finite N 
effects can be described by a multiplicative correction. 

Note that the upper edge b(s) = (s+l) s+1 /s s of the FC 
distribution tt( s > (x) for large matrices determines the size 
of the largest eigenvalue X m ax of the reduced density ma- 
trix p of size N. In the case of the structureless ensemble 
of random pure states on N x N system, corresponding to 
s = 1, one has 6(1) = 4 so that A mol « 4/N. For states 
obtained by projection of a 2s partite system on maxi- 
mally entangled states, as described in the previous sec- 
tion, the largest component behaves as X ma x ~ b(s)/N. 

The number A maa; , equal to the largest component of 
the Schmidt vector of a random pure state on the bi- 
partite system, can be used to measure the degree of 
quantum entanglement. For instance, for a bi-partite 
system, the 'geometric measure' of entanglement, related 
to the distance to the closest separable state (in sense of 
the natural Fubini-Study distance) reads 43], E g {\<j>)) — 
— In A max . This quantity can also be considered as the 
Chebyshev entropy Soo - the generalized Renyi entropy 



S q = lnTrp 9 in the limit q — > oo 0. 

Thus the right edge b(s) of the support of the spec- 
tral density for the reduced state p = Tr b\4>) deter- 
mines the geometric measure of entanglement of the cor- 
responding random pure states. In the case of struc- 
tureless random pure states, related to the Marchenko- 
Pastur distribution one becomes an asymptotic expres- 
sion (E g (\(p)) = — ln(4/iV). In the general case of random 
state corresponding the the FC distribution n^(x) the 
typical value of the entropy reads 



(5oo) s =kxN-kxb(s) = InJV + sins - (s + 1) ln(s + 1). 

( 41 ) 

The larger value of s, the smaller the Chebyshev entropy 
Soo, and the less entangled a typical random state ob- 
tained by the projection of the 2s partite system. 

The average von Neumann entropy, Si = — Trplnp, 
of mixed states of size N generated according to the FC 
distribution reads (S(p)) s = lnN— J^i 2 J @- ^ ne sec " 
ond moment of the FC distribution given in (|35l) implies 
the asymptotic average purity (Trp 2 ) s rj (s+ l)/N - the 
larger number s, the more pure the typical mixed state 
generated by a projection onto s — 1 maximally entangled 
states. 
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VII. CONCLUDING REMARKS 

In this work we analyzed structured ensembles of ran- 
dom pure states on composite systems. They are defined 
with respect to a given decomposition of the entire sys- 
tem into its subsystems, what induces a concrete tensor 
product structure in the Hilbcrt space. The structured 
ensembles are thus invariant with respect to the group of 
local unitary transformations. 

Performing a partial trace over selected subsystems one 
obtains an ensemble of random mixed states defined on 
the remaining subsystems. The particular ensemble de- 
pends thus on the number of systems traced out and on 
the way the initial random pure states are prepared. 

Quantum states obtained by the partial trace of a su- 
perposition of k maximally entangled pure states of the 
bi-partite system involve the sum of k random unitary 
matrices. To generate states which involve a product 
of an arbitrary number of s matrices one needs to con- 
sider a system consisting of 2s subsystems, in which an 
orthogonal measurement is performed in the product of 
s — 1 maximally entangled bases. Selected ensembles of 
random states, recipe to generate numerically the cor- 
responding density matrices and some properties of the 
distribution of the Schmidt coefficients are collected in 
Table 1. 

We are going to conclude this work by writing down 
a more general class of structured random states, which 



contains all particular cases discussed in the paper and 
listed in Table 1. Consider the following ensemble of non- 
hermitian random matrices parametrized by an arbitrary 
fc-dimensional probability vector p — {p\, . . . ,pk\ and a 
non-negative integer s, 



W kii 



[p\Ui +p 2 U 2 



in.-i;>, (;,■■■(;.. (42) 



Here U\, . . .U k denote k independent random unitary 
matrix distributed according to the Haar measure on 
U(N), while Gi,...G s are independent square random 
matrices of size N from the complex Ginibre ensem- 
ble. Random density matrix is obtained as a normalized 
Wishart-like matrix, 



W k .W, 



PK 



Tr(W k>s Wl s ) 



(43) 



Note that any particular ensemble from the above class 
can be physically realized by taking a superposition of k 
random pure states weighted by the vector p. Each pure 
state is defined on the system containing 2s subsystems. 
Performing a measurement in the product of (s — 1) max- 
imally entangled bases one gets a random pure state of 
the desired structured ensemble. Eventually, averaging 
over the last subsystem one arrives at the mixed state 



k 


s 


matrix W 


distribution P(x) 


singularity at x — > 


support [a, b] 


M 2 


mean entropy 


1 







8(1) = 7T<°> 




{1} 


1 





2 





Ux+U 2 


arcsine 




[0,2] 


3/2 


In 2- 1 « -0.307 


3 





Ui+U 2 + U 3 


3 entangled states 




[0,2f] 


5/3 


» -0.378 


4 





Ul + U 2 + U 3 + Ua 


4 entangled states 




[0,3] 


7/8 


« -0.411 


1 


1 


G ~ UG 


Marchenko-Pastur 7P 1 ' 


x- 1 ' 2 


[0,4] 


2 


-1/2 = -0.5 


2 


1 


(t/i + U 2 )G 


Bures 




[0, 3V3] 


5/2 


- In 2 « -0.693 


1 


2 


G\G 2 


Fuss-Catalan 7r' 2 ' 




[0,6f] 


3 


-5/6 w -0.833 


1 
1 


s 


Gi • ■ • G s 


Fuss-Catalan 7r' s ' 


x~ s/(s+1) 


[o,( S +rr +i /s s ] 


8+1 


\ps+l 1 



TABLE I: Reduction of pure states from structured ensembles leads to random mixed states of the form p — WW^ /TrWW^ . 
Random matrix W is constructed out of random unitary matrices Ui distributed according to the Haar measure and/or 
(independent) random Ginibre matrices Gj of a given size N. Asymptotic distribution P(x) of the density of a rescaled 
eigenvalue x = iVA of p for N — > 00 is characterized by the singularity at 0, its support [a,b], the second moment M 2 
determining the average purity (Trp 2 ) = M 2 /N and the mean entropy, — x In xP(x)dx, according to which the table is 
ordered. 



Consider first the case of a uniform probability vector, 
Pi = 1/k for i = 1, . . . , k. For k — 1 one obtains ensem- 
bles leading to the Fuss-Catalan distributions tt^ s \ which 
in the case s = 1 reduces to the Marchenko-Pastur distri- 
bution. Taking k — 2 and s = one obtains the arcsine 
ensemble ([7]), while for larger k one obtains the distribu- 



tions (|10|) . which converge to ttW in the limit k —¥ 00. 
Moreover, the case k — 2 and s = 1 corresponds to the 
Bures ensemble (|23|). Thus the case k = 2 and arbitrary 
s can be called higher order Bures ensemble. 

In a more general case, taking an arbitrary probability 
vector p and varying the weights in a continuous man- 
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ner one can study transition between given structured 
ensembles. For instance, by fixing the parameter s, set- 
ting k = 2 and varying the weight p2 = 1 — pi one de- 
fines a continuous interpolation between the higher or- 
der Bures ensemble and the Fuss-Catalan ensemble. We 
have shown therefore that having in our disposal sim- 
ple algorithms to generate random unitary and random 
Ginibre matrices we can construct a wide class of ensem- 
bles of random quantum states. Furthermore, we provide 
constructive physical recipe to generate such states by 
means of generic two-particle interaction, superposition 
of states, selective measurements in maximally entangled 
basis and performing averages over certain subsystems. 

As discussed in Appendix B it is also possible to intro- 
duce analogous ensembles of real random density matri- 
ces. Physically this corresponds to imposing restrictions 
on the class of the interactions used to generate random 
pure states. In contrast with the complex case, the en- 
semble based on square real Ginibre matrices does not 
lead to the Bures measure in the space of real states. 
To achieve such a measure one needs to generalize the 
ensemble even further to allow also rectangular Ginibre 
matrices 15j. In physical terms this implies that the di- 



mension of the principal system and the auxiliary system 
have to be different in this case. 

The notion of random quantum states is closely re- 
lated with the concept of random quantum maps. Due 
to the Jamiolkowski isomorphism any quantum operation 
$ acting on density matrices of size N can be represented 
by a state on the extended Hilbert space [S] , 



(44) 



Here \ip + ) = \y^^2jLi \ denotes the maximally en- 
tangled state from the bipartite Hilbert space H = 
%n ® Hn- Any state a on the composed Hilbert space 
% = %a ®T~Lb defines a completely positive, trace pre- 
serving map provided the following partial trace condi- 
tion is satisfied, Tr ao{$) = 1/N. 

It is possible to impose this partial trace condition on 
an arbitrary state u> acting on %. To this end one finds 
the reduced state Y := T"aw which is positive and al- 
lows one to take its its s quar e root y/Y, and writes the 
normalized cognate state [44j |. 



a=—(l® -=)uj(l ® -7=1 ; 
N\ Jy) \ Jy J 



(45) 



The required property, Tr^tr = 1/N, is satisfied by con- 
struction, so the state a represents a quantum operation. 
As the matrix elements of the corresponding superopera- 
tor $ can readily be obtained by reshuffling Q the entries 
of the density matrix a, any random state uj determines 
by (|45|) and (|44]) a quantum operation. Therefore any 
ensemble of random states introduced in this paper, ap- 
plied for bi-partite, N x N systems determines the corre- 
sponding ensemble of random operations. For instance, 
the induced measure with K — N 2 corresponds to the 
flat measure in the space of quantum operations [44], EH , 



but other ensembles of random states can be also applied 
to generate random quantum operations [9] . 

The present study on ensembles of random states 
should be concluded with a remark that apart of the 
methods developed in this paper several other approaches 
are advocated in the literature. In very recent papers 
[46l |47| the authors follow a statistical approach intro- 
ducing a partition function which leads to a generaliza- 
tion of the Hilbert-Schmidt measure. Varying the pa- 
rameter of the model, which corresponds to the inverse 
temperature, they demonstrate a phase transition during 
an interpolation between Marchenko-Pastur and semicir- 
cle distribution of spectral density. In another recent 
approach Garnerone et al. study statistical properties of 
random matrix product states [48j . which are obtained 
out of products of truncated random unitary matrices. 
Although these models of random states do differ from 
the one presented in this work, a possible links and re- 
lations between results obtained in these approaches is 
currently under investigation. 
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Appendix A: Sum of k random unitaries and the 
distribution vu 

We compute, in the limit of large matrix size N — > 
oo, the asymptotic eigenvalue distribution of the random 
matrix 



v k = ^(U 1 + U 2 + 



+ u k )(u* + u; + -- ■ + u^) > (Ai) 



where U\ , . . . , U k are N x N random independent Haar 
unitary matrices. Obviously, this is equivalent to com- 
puting the singular value distribution of A: -1 / 2 Yli=i Ui- 
For now, we forget about the normalization pre-factor 
and we put W k = kV k . It is a well known result in free 
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probability theory that independent large unitary matri- 
ces are free from each other (and, very importantly but 
not of interest here, they are also free from deterministic 
diagonal matrices): 

Theorem A.l ([49fl and [Hoj]) Let U 1>N , . . . , U k . N £ 
U(N) be k independent Haar unitary random matrices. 
Then, as N — > oo, 

Ui, N ,...,U k ,N * dlStr > ui, . . . ,u k , (A2) 

where u\ , . . . , u k are free Haar unitary elements in a non- 
commutative W* -probability space (A4,t). 



In [52j , it is shown that the Brown measure of w is a ro- 
tationally invariant measure, supported on the centered 
disk of radius y/k with radial density 

M r ) = ^i fc ~2?2 ' 0<r<Vk. (A8) 
TT{k z — r z ) z 

With the proper fc -1 / 2 rescaling, it is easy to see that the 
above Brown measure converges to the uniform measure 
on the unit disk. Hence, we recover the Ginibre behav- 
ior in the limit k — s- 00 (one has to take first the limit 
N — > 00). Let us add that a more general study on sta- 
tistical properties of a sum of random unitary matrices 
was recently presented by Jarosz [531 ] . 



Hence, computing the limit distribution of W k 
amounts to understanding the distribution of a sum of 
k free Haar unitary elements in a von Neumann alge- 
bra w k = (u\ + • • ■ + Uk)(u* + ■ ■ ■ + u^). This problem 
has been related to random walks on fc-regular trees by 
Kesten [5l|. Indeed, the number of alternating words of 
length 2p in the letters Ui,u* which reduce to the unit is 
bijectively the same as the number of walks of length 2p 
on the /c-regular tree beginning and ending at some fixed 
vertex. Using standard formulas for the number of such 
walks, we can deduce moment information for k = 2 



lim E 

JV-s-oo 



N 2 



(A3) 



and a moment generating function in the general case 



2(fc- 1) 



p=0 



2 + ky/l-4(k-l)z 



(A4) 



From the last formula, using Cauchy transform tech- 
niques, one can easily deduce the probability density 
function of the distribution of w: 



d/j k (x) 



k \J^{k - Tjx - x 2 



[o,4(k-i)]{x)dx. (A5) 



This result has been obtained by Haagerup and Larsen 
in [52j, Example 5.3. The authors were interested in the 
Brown measure of the non-normal element w = u± +• • • + 
u k . They found that the distribution of \w\ is given by 



. , k JA{k-l)-x 2 



2Vk^T]( x ) dx - ( A6 ) 



One can easily recover equation (|A5[) by using w = ww* 
and by noticing that fi k = A*fc sq, where sq is the square 
function sq(a;) = x 2 . 

If v k is the distribution of the rescaled element v = 
w/k, then one arrives at the desired distribution function, 



dv k (x 



1 y/4k(k- l)x-k 2 x 2 

) — ~ ~ ~ l[0,4fci]( 



2n 



kx 



'x)dx. (A7) 



Appendix B: Real random states, real Ginibre and 
random orthogonal matrices 

Although most often one considers complex density 
matrices, it is also interesting to study quantum states 
described by real density matrices. The dimensionality 
of the set of real states on Hn is N(N + l)/2 — 1, so 
its geometry is easier to study than that of the N 2 — 1 
dimensional set of complex states Q. For instance, the 
set of real states of a qubit forms a two dimensional disk, 
which can be considered as a cross-section of the three 
dimensional Bloch ball of complex states. Euclidean vol- 
ume of the set of real density matrices of size N was 
derived in [54 |. while the corresponding measure can be 
derived from the real Ginibre ensemble. 

In this appendix we define ensembles of real states 
based on random orthogonal matrices and real Ginibre 
ensemble and show that the level density does not differ 
from the complex case. To this end we formulate two 
lemmas. 

Lemma B.l Consider k independent orthogonal matri- 
ces 0\ , . . . O k , distributed according to the Haar measure 
on O(N) and define a normalized real density matrix 



Port — 



{Ox 



O k )(Of 



+ ol) 



Tr(Ox 



O k ){0{ 



ol) 



(Bl) 



Then for large N its spectral density is described by the 
distribution v k given in ilO\) . 



This lemma follows directly from the fact that the mo- 
ments of orthogonal and unitary random matrices have 
the same behavior for large matrix size N, since random 
orthogonal matrices are asymptotically free [55| . 

Lemma B.2 Consider s independent random matrices 
R\ , . . . R s taken from the real Ginibre ensemble of square 
matrices of size N . Define a normalized real density ma- 
trix 

R1R2 ■ ■ ■ R s (RiR2 ■ ■ ■ Rs) T 



PR 



Tv[R l R 2 ---R s {R l R 2 ---R s ) T ] 



(B2) 



Then for large N its spectral density is described by the 
Fuss- Catalan distribution tt^ given in 
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To prove this one needs to show that in the case of 
large matrix size N the moments of this distribution are 
indeed given by the Fuss-Catalan numbers (1551 , exactly 
as for the product of complex Ginibre matrices. This 
follows from the interpretation of the Wick formula for 
random matrices in terms of maps gluing and from the 
fact that leading terms have to be non-crossing, therefore 
orientable as for complex Gaussian matrices - see e.g. 

mm. 

It is natural to combine both definitions and defining a 



more general ensemble of real density matrices, each ob- 
tained out of k random orthogonal matrices and s square 
real Ginibre matrices in a direct analogy to eq. (|42[) . 

Let us close this section with a remark that the dif- 
ferences between the real and complex case can be sig- 
nificant in some cases. For instance, to get a real state 
distributed according to the Bures measure one needs to 
use a symmetric random unitary matrix and a rectangu- 
lar, N x (N + 1), real Ginibre matrix, while the complex 
Bures state is obtained from a random unitar y a nd a 
square matrix of the complex Ginibre ensemble [15| . 
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